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SUMMARY 


PROBLEM 


The  increased  emphasis  on  the  use  of  simulators  for  train¬ 
ing  of  infrared  (IR)  and  low  light  level  television  (LLLTV) 
systems  has  presented  some  technical  simulation  problems. 
One  such  problem  is  target  identification  which  requires 
the  need  for  proper  textural  cues,  especially  on  low-level, 
long-range  flight  profiles.  These  textural  cues  aid  in 
providing  Increased  simulator  fidelity  through  enhanced 
realism  and  better  discrimination  between  target  types. 


APPROACH 


Computer  image  generation  (CIG)  techniques  offer  solutions 
to  these  simulation  problems.  However,  the  capacity  of 
typical  edge-type  algorithms  for  textural  detail  rapidly 
outstrips  system  capabilities  for  real-time  operations. 

One  alternate  technique  is  the  use  of  cyclic  codes  for  pro¬ 
viding  surface  detail  that  is  spatially  correctable  and 
provides  different  levels  of  detail  as  a  function  of  range. 
This  technique  was  used  on  a  previous  commercial  contract 
to  generate  a  complete  seascape  in  real  time  on  a  TV  raster- 
type  display.  For  this  AFHRL  study,  hardware  was  modeled  in 
software  to  provide  a  non-real-time  simulation  that  was 
applicable  to  both  real-time  operations  and  present  CIG- 
type  simulators.  CIG-type  coordinate  trans formati ons , 
perspective  equations,  clipping,  and  edge-smoothing  algo¬ 
rithms  were  written  for  operation  of  the  simulation  on 
the  AFHRL  Advanced  CIG  Techniques  Evaluation  System  (ACTES) 
facility  at  Wri ght-Patterson  AFB. 

A  data  base  definition  schema  was  developed  using  the  com¬ 
pact  notation  of  cyclic  codes,  their  theoretical  proper¬ 
ties  and  algorithms  for  varying  the  code  as  a  function  of 
geographic  position.  The  schema  is  based  upon  non-edge- 
type  algorithms. 


RESULTS 

A  simulation  program  was  written  and  made  operational 
on  contractor  faciliites  and  the  Xerox  Sigma  5  and  ACTES 
facilities  at  WPAFB.  The  simulation  was  displayed  on  the 
ACTES  display  using  Its  offline  data  storage  and  display 
capabilities.  The  simulation  demonstrated  the  types  of 
texture  that  could  be  generated  from  a  cyclic  code  data 
base.  Variation  In  levels  of  detail  from  background  to 
foreground, as  well  as  spatial  coherence,  was  demonstrated 


by  running  various  scenarios.  Geometric  and  texturing 
aspects  were  demonstrated  for  various  fields  of  view  and 
aircraft  positions  with  respect  to  the  textured  terrain 
area*  showing  the  applicability  of  the  techniques  to  both 
sensor  and  visual  type  displays. 


CONCLUSIONS 

The  implementation  on  the  ACTES  demonstrates  that  cyclic 
coding  techniques  are  feasible  from  mathematical  and  soft¬ 
ware  programming  viewpoints  for  simulators  typifying  oper¬ 
ational  IR  and  LLLTV  systems.  Texture  realism  was  not 
validated  as  it  was  not  a  part  of  this  study.  Future 
efforts  should  include  three-dimensional  aircraft  motion, 
texture  realism  and  validation  of  the  real-time  aspects  of 
the  implemented  software  models  and  algorithms. 


PREFACE 
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Air  Force  Human  Resources  Laboratory,  Wri ght-Patterson 
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Digital  Communications,  214  Flynn  Avenue,  Moorestown,  New 
Jersey  08057.  Mr.  Joseph  F.  deSpautz  was  the  principal 
Investigator  for  Applied  Digital  Communications.  Captain 
Michael  L.  Ingalls  was  the  contract  manager  for  the  Air 
Force  Human  Resources  Laboratory.  Mr.  Michael  Nicol  of 
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conversion  to  the  Sigma  5  computer  and  in  running  the  simu¬ 
lation. 
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1.0  INTRODUCTION 

The  purpose  of  this  study  was  to  demonstrate  the  feasibility 
of  using  pseudo  random  noise  (PRN)  codes  for  surface  tex¬ 
turing  in  flight  training  simulators.  To  accomplish  this, 
a  software  simulation  was  developed  to  produce  static  and 
dynamic  views  of  a  flat  surface  from  varying  viewpoints  in 
simulated  flight,  illustrating  changes  in  texture  detail 
with  changes  in  perspective.  Included  among  the  system 
requirements  were  the  ability  to  control  amplitude  and 
spacing  of  tonal  variations  and  to  demonstrate  flexibility 
in  generating  a  variety  of  texture  patterns. 

Previous  applications,  using  cyclic  codes  for  image  enhance¬ 
ment,  have  used  the  randomness  or  incoherent  code  properties 
for  white  noise  generation  that  is  added  to  the  scene.  This 
type  of  texturing  is  not  correlated  TV  frame  to  TV  frame, 
and,  therefore,  is  not  usable  for  simulation  realism  or 
training  cueing.  This  study  approach  associates  code  states 
to  geographic  position  (or  entire  scene)  and  uses  the  code 
properties  and  methods  for  efficiently  generating  them, 
thereby  allowing  line-to-line  and  TV  frame- to- frame  corre¬ 
lation. 

The  texture  is  calculated  in  a  deterministic  manner  and 
overlays  a  constant  tone  ground  area.  Scene  fidelity  will 
be  demonstrated  by  showing  the  variation  in  textural  de¬ 
tail  as  a  function  of  sensor  altitude  and  position. 

The  implementation  of  PRN  coding  techniques  has  already 
been  shown  in  company  programs,  including  a  real-time  sea¬ 
scape  generator.  Cyclic  code  generators  can  operate  at  TV 
raster  generation  speeds.  PRN  structures,  tonal  variation, 
refresh  rate,  etc.  can  all  be  controlled  in  real  time  with 
appropriate  hardware. 

Large  area  texturing  can  be  readily  constructed  with  a 
minimum  of  storage  requirements.  The  initial  texture 
definitions  need  only  be  kept  in  memory  for  each  of  the 
PRN  generators,  and  these  values  can  be  eliminated  with 
the  construction  of  a  hardware  function  to  calculate  the 
individual  starting  states  in  real  time. 

To  demonstrate  the  texturing  capabilities  of  using  PRN 
techniques  for  surface  texturing,  scenarios  of  level  flight 
over  a  rectangular  terrain  area  that  has  a  texture  were 
recorded  for  playback  on  the  AFHRL  ACTES  facility.  To 
permit  a  variety  of  test  cases,  parameters  including 
terrain  size  and  location,  PRN  structure,  aircraft  altitude 
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and  velocity,  and  field  of  view  are  under  interactive 
control  of  the  user.  Incorporated  in  the  simulation  are 
the  typical  operations  of  spatial-image  plane  transforma¬ 
tions,  field-of-view  clipping,  perspective  transformations, 
pixel  intensity  calculations,  edge  smoothing,  and  display 
formatting. 

The  simulation  demonstrates  several  types  of  textures  which 
are  possible  with  the  cyclic  code  data  base.  Control  of 
textures  is  achieved  in  two  ways:  (a)  the  use  of  repeating 
patterns  that  are  generated  by  cyclic  codes  which  are  called 
subterrains  permits  reinitialization  of  the  basic  PRN  pat¬ 
tern  at  various  points  to  introduce  desired  discontinuities, 
and  (b)  the  degree  of  texture  detail  can  be  controlled  by 
reducing  the  range  of  possible  gray  levels  to  any  desired 
number. 

Geometric  and  texturing  aspects  have  been  demonstrated  from 
various  fields  of  view  of  the  sensor  display  and  aircraft 
positions  relative  to  the  textured  surface.  As  the  range 
of  the  aircraft  to  the  data  base  is  decreased,  the  expected 
increase  in  resolvable  detail  is  demonstrated  quite  accur¬ 
ately.  Alterations  in  display  field  of  view  show  that  the 
technique  is  potentially  applicable  to  visual  as  well  as 
sensor  simulations. 


2.0  TECHNICAL  APPROACH 


A  number  of  different  techniques  have  been  investigated  to 
provide  surface  detail  or  terrain  texturina  in  computer 
image  generation  (CIG)  for  simulation  of  operational  infra¬ 
red  (IR)  and  low  light  level  television  (LLLTV)  systems. 
These  techniques  include  cyclic  coding  algorithms  for  tonal 
generation. 

Past  company  programs  have  demonstrated  that  compact  cyclic 
codes  could  be  used  to  provide  a  coherent  real-time  surface 
texturing  of  dynamic  seascape  scenes.  This  is  contrasted 
to  other  applications  of  cyclic  code  techniques  that  provide 
texturing  using  the  random  or  incoherent  code  properties 
(i.e.,  white  noise  generator).  This  latter  type  of  tex¬ 
turing  is  not  spatially  correlated  and,  therefore,  is  not 
valid  for  display  realism. 

The  approach  used  in  this  study  was 

1.  To  associate  code  states  to  geographic  position 

to  provide  spatial  correlation  or  visual  coherence. 

2.  To  develop  software  models  of  rtal-time  hardware 
techniques  for  the  generation  of  texture  for 
presentation  on  TV  raster-type  displays. 

3.  To  develop  a  data  base  definition  without  edges 
for  variable  surface  texture  that  is  parameter 
controlled. 

4.  To  develop  a  non-real-time  software  simulation 
for  execution  in  the  AFHRL  ACTES  facility. 

5.  To  demonstrate  various  texturing  concepts  by 
running  and  recording  scenarios  on  existing 
AFHRL  facilities. 

The  system  flow  diagram  of  the  simulation  developed  is 
shown  in  Figure  1  and  includes: 

1.  Data  Base  Generation. 

2.  Data  Input/Output  and  Simulation  Control. 

3.  Aircraft  Sensor  Positioning. 

4.  Spatial  Transformations. 

5.  Perspective  Transformations. 


I  I 
| DATA  BASE  , 
DEFINITION 


PEN 

AVERAGING 


{ FROM \ 

(controller  ) 


DATJ 
AND  SI1 
C01 

\  I/O 

'lULATOR 

4TROL 

aircraft 

POSITIONING 
AND  SPATIAL 
TRANSFORMATION 

TV  PIXEL  TONE 
CALCULATION 

1 

VIDEO 

ASSEI 

LINE 

•ELY 

BASIC 

SIMULATION 

PARAMETERS 


Figure  1.  System  Flow  Diagram 
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6.  TV  Raster  Pixel  Tone  Calculation. 


7.  Video  Line  Assembly. 

Described  in  this  section  are  the  technical  approach  and 
mathematical  techniques  used  in  the  simulation  to  develop 
the  dynamic  TV  frame  data  which  was  stored  on  magnetic 
tape.  The  data  tape  uses  the  offline  ACTES  facility  for 
visual  display. 

A  discussion  of  data  base  construction  is  presented.  Since 
the  terrain  texture  is  defined  by  using  PRN  codes  to  repre¬ 
sent  texture  tones,  an  analysis  of  PRN  codes  and  their 
generation  is  included.  The  mapping  of  these  codes  onto 
a  terrain  surface  is  described  which  defines  the  terrain 
structure. 

The  simulation  includes  the  coordinate  transformations  for 
independent  aircraft  motion,  clipping  of  terrain  for  a 
variable  field  of  view  (POV)  screen  with  selectable  para¬ 
meters  for  image  display  characteristics,  and  algorithms 
for  determining  the  visible  portion  of  the  terrain  within 
each  TV  pixel.  The  initial  study  was  limited  to  horizon¬ 
tal  flybys  over  the  terrain,  and  rotational  effects  were 
investigated  under  a  contract  amendment. 

An  algorithm  associates  the  correct  area  of  the  terrain 
with  a  given  pixel  by  projecting  the  pixel  definition  onto 
the  terrain  and  determining  the  location  and  size  of  the 
projection.  The  resulting  composite  area  tone  is  the 
image  plane  value  that  is  to  be  displayed  if  the  terrain 
is  not  blocked  by  any  other  objects. 

To  provide  different  levels  of  texture  detail,  multiple 
PRN  values  fall  within  a  pixel;  thus,  a  technique  is  pro¬ 
vided  for  averaging  these  values  with  a  minimum  of  calcu¬ 
lation.  For  cases  in  which  a  terrain  edge  or  boundary 
falls  within  a  pixel,  edge  smoothing  techniques  are  applied. 

The  simulation  was  optimized  in  terms  of  execution  time  and 
amount  of  texturing  control.  Integer  mathematics  was  im¬ 
plemented  to  eliminate  floating  point  arithmetic  to  improve 
run  times.  Control  of  textural  tone  attenuation  was  made 
a  program  control  parameter. 
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2.1  DATA  BASE  DEFINITION 

2.1.1  CYCLIC  CODE  CONSIDERATIONS 

One  of  the  most  important  advances  in  error  correcting 
code  development  has  been  concerned  with  cyclic  codes.  A 
cyclic  code  is  one  in  which  each  element  of  a  set  is  ob¬ 
tained  by  cyclically  shifting  each  predecessor  one  unit. 

Cyclic  codes  are  of  considerable  importance  in  CIG  enhance¬ 
ments  due  to  their 

1.  Theoretical  properties. 

2.  Ease  of  generation. 

3.  Ease  of  control. 

The  cyclic  codes  chosen  for  investigation  in  this  study 
are  PRN  codes  having  the  following  properties: 

1.  Maximal  length 

2 .  Polynomial  degree 

3.  Primitive  polynomial  characteristics 

4 .  Compact  polynomial  notation 

Maximal  length  refers  to  the  code  set  having  maximum  peri¬ 
odicity;  that  is,  the  number  of  unique  states  which  the 
code  generates  before  the  code  is  repeated. 

Polynomial  degree  refers  to  total  length  of  the  code  se¬ 
quence  and  is  an  indication  of  the  amount  of  granularity 
or  variation  in  a  sequence. 

A  primitive  polynomial  indicates  a  unique  code.  A  primi¬ 
tive  code  is  not  reducible  to  two  other  codes. 

Polynomials  offer  the  advantage  of  compact  notation.  For 
a  maximal  length,  primitive  polynomial  of  degree  m,  2™-! 
code  states  are  generated. 

2.1.2  LINEAR  AND  NON-LINEAR  FUNCTIONS 

Shift  registers  with  linear  feedback  functions  are  the 
most  amenable  to  analysis.  Linear  Boolean  functions  are 
analogous  to  linear  algebraic  functions.  Any  linear  al¬ 
gebraic  equation 

ax  +  by  =  c 

may  be  solved  for  y  once  x  is  given. 
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The  Boolean  equation 
x  +  y  =  1 

is  a  non-linear  function  of  x  and  y.  The  +  sign  here  de¬ 
notes  logical  addition  (OR) .  Given  a  state  value  of 

x  =  1 

y  is  indeterminate.  Since  y  cannot  be  determined  for  both 
state  values  of  x,  the  binary  operation  of  logical  addition 
is  referred  to  as  information  losing.  A  linear  operation 
must  be  information  preserving. 

Logical  multiplication  is  also  non-linear.  Given  the 
Boolean  equation 

x  •  y  =  0 

and  the  state  value 
x  =  0 

y  is  indeterminate. 

The  Boolean  operation  of  addition  modulo-2  (exclusive  OR) 
is  linear  and  information  preserving.  Given  the  Boolean 
equation 

x  ®  y  =  1 

or  x  ©  y  =  0 

y  may  be  determined  once  x  is  given. 

Only  linear  functions  will  be  considered  here. 

2.1.3  MATRIX  DESCRIPTION  OF  CYCLIC  CODES 

Consider  the  generalized  r-stage  feedback  shift  register 
with  a  linear  feedback  function  (Figure  2) .  The  stages  of 
the  register  are  designated  Sg,  Si,  S2,  ...»  S^..  The  out¬ 
put  of  stage  S^  is  connected  to  the  input  of  tne  modulo-2 
summer  when  C<  *  1.  Cg  is  always  one,  otherwise  fewer  than 
r  stages  would  be  in  use.  The  external  input,  e,  is  a 
Boolean  constant. 

Let  x^  represent  the  current  state  and  X^  the  next  state 
of  stage  S^.  The  next  state  X^  of  each  stage  may  be  ex¬ 
pressed  as  a  linear  Boolean  function  of  the  present  state 
x^  of  one  or  more  stages. 
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Generalized  r-Stage  Feedback  Shift  Register 


Xr  *  crxr  ©  cr_1  xr_1  ©  . 


.  9  Cj^Xj^  e  c0x0  e  e 


Vi  =  xr 
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This  system  of  equations  may  be  expressed  as: 
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or  X  =  Tx  9  L 


The  r  x  r  Boolean  matrix,  T,  is  non-singular  since  its 
row  (column)  vectors  are  linearly  independent.  T  repre¬ 
sents  the  linear  transformation  of  an  r  component  vector 
(present  state)  into  another  r  component  vector  (next 
state).  L  represents  a . translation.  When  e  •  1,  the 
modulo- 2  sum  of  the  vector  L  with  Tx  represents  the  com¬ 
plement  of  the  bit  fed  back. 

A  linear  transformation  T  followed  by  a  translation  L  is 
called  an  affine  transformation.  Every  translation  is 
one-to-one  and  has  an  inverse.  The  linear  transformation 
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T  is  one-to-one  and  has  am  inverse.  Hence,  the  affine 
transformation  Tx  ®  L  is  one-to-one  and  has  an  inverse. 
Because  of  this  property,  each  state  has  a  unique  successor 

Of  primary  interest  are  those  feedback  combinations  which 
yield  maximal  length  sequences.  Those  maximal-length 
cyclic  codes  with  e  =  0  have  been  investigated. 

For  the  case  e  =  0, 

X  =  Tx 

An  initial  state  x  will  be  followed  by  the  states 
Tv  T^x  T*3  •  •  •  s  Tk+1  —  ip 

The  smallest  value  of  k  for  which  T^  =  I  is  the  length  of 
the  longest  possible  cycle.  In  general,  for  maximal 
length  sequences 

k  -  2m  -  1 

where  m  is  the  degree  of  the  feedback  polynomial. 

2.1.4  APPLICATION  TO  DATA  BASE  DEFINITION 

The  data  base  generation  algorithm  to  create  a  textured 
grid  steps  the  initial  starting  state  p^  using  the  appro¬ 
priate  code  generator  to  yield  texture  values  Pi+^,  Pj.+2' 
etc.  The  total  textured  area  of  the  terrain  is  generated 
by  altering  the  starting  state  per  grid  line  and  perform¬ 
ing  the  code  generation  algorithm  with  the  new  starting 
state.  An  example  is  to  cause  the  starting  state  to  pre- 
cess  by  one  state  for  each  line,  thereby  producing  a  skewed 
field.  This  precession  by  one  code  state,  yielding  a 
symmetric  configuration  of  codes,  has  a  number  of  mathe¬ 
matical  properties  which  are  presented  later  in  this 
section.  It  is  the  data  base  texture  procedure  used  in 
this  investigation. 

The  software  which  generates  these  cyclic  codes  was  de¬ 
signed  to  simulate  the  actual  hardware  algorithm.  A  pseudo 
shift  register  was  implemented  using  a  table  lookup  algo¬ 
rithm  to  decrease  software  execution  time. 

2.1.5  TERRAIN  STRUCTURE 

The  size  and  shape  of  the  terrain  that  is  contained  in 
each  pixel  projection  determines  the  PRN  tones  to  be  com¬ 
bined  or  averaged  to  form  the  pixel  tone.  There  are 
potentially  three  types  of  combinational  operations  that 
have  to  be  considered  to  provide  proper  pixel  tone.  These 
operations  are  a  function  of  data  base  quantization  and 
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object-image  plane  geometry  and  include  the  following 
cases: 

1.  One  PRN  code  overlays  more  than  one  pixel. 

2.  One  PRN  occupies  one  pixel. 

3.  Many  PRN  codes  and  portions  thereof  occupy 
a  pixel. 

The  latter  case  is  the  principal  case  investigated  in  this 
study  as  it  shows  the  texture  variation  as  a  function  of 
range  and  is  representative  of  realistic  conditions. 

Figure  3  shows  the  representative  geometry  of  the  terrain 
structure.  The  data  base  is  a  rectangular  area  identified 
by  its  vertices  in  the  X,Y  plane.  The  vertices  form  edges 
for  use  in  the  perspective  and  clipping  algorithms.  The 
data  base  is  quantized  into  a  grid  in  which  each  grid  ele¬ 
ment  defines  a  "minimum  texture  area"  (MTA) .  Each  of  these 
MTA  elements  is  associated  with  a  PRN  value  representing 
the  texture  tone.  The  terrain  data  base  lying  on  the  sur¬ 
rounding  ground  has  a  higher  priority  than  the  surrounding 
ground  and,  therefore,  blocks  it  in  the  image  plane.  In 
the  simulator,  the  data  base  size  and  MTA  size  are  simu¬ 
lation  variables  controlled  by  input  parameters. 

Each  PRN  value  is  associated  with  a  tonal  intensity  for 
each  MTA.  A  row  of  MTAs  parallel  to  the  Y-axis  forms  a 
terrain  line.  The  assignment  of  PRN  values  to  a  terrain 
line  is  made  by  specifying  a  "starting  state”  for  the  first 
MTA  on  that  line.  A  starting  state  is  an  index  into  the 
sequence  of  PRN  values;  it  is  used  to  associate  a  particu¬ 
lar  PRN  value  with  the  MTA  to  which  it  is  assigned. 

Within  a  group  of  adjacent  terrain  lines  called  a  "sub¬ 
terrain,"  starting  states  are  increased  by  one  from  line 
to  line.  This  arrangement  permits  the  use  of  techniques 
for  estimating  the  average  tone  of  multiple  MTAs.  By 
varying  the  size  of  subterrains,  as  well  as  the  starting 
state  differences  between  adjacent  subterrains,  the  over¬ 
all  texture  of  the  entire  terrain  can  be  altered. 

The  subterrain  definition  is  illustrated  in  Figure  3.  The 
first  subterrain  has  a  starting  state  of  1  assigned  to  the 
first  terrain  line.  Within  it,  starting  states  differ  by 
a  slip  rate  of  one.  The  final  starting  state  of  the  first 
subterrain  is  five.  For  the  next  subterrain,  the  initial 
starting  state  was  slipped  5  states.  Within  the  sub¬ 
terrain,  however,  the  state  of  each  terrain  line  again  is 
incremented  by  one  state  continuing  the  same  pattern  as 
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■•-Terrain  Line-*- 


Data  Base  Definition 


in  the  first  subterrain  but  with  different  codes. 

Using  this  feature,  different  patterns  could  be  generated 
from  the  same  PRN  code  and  simulation  algorithms,  thereby 
showing  the  capabilities  of  the  texturing  concepts. 

The  actual  size  and  shape  of  the  texture  area  for  the 
straight-and- level  flight  scenarios  that  were  demonstrated 
in  this  study  are  long  and  narrow  quadrilaterals.  Because 
the  projections  of  the  pixels  that  form  a  TV  line  are  para¬ 
llel  to  the  fixed  coordinate  system,  the  actual  texture 
area  is  approximated  by  rectangles  whose  coordinates  are 
determined  from  the  perspective  equations. 

The  approximation  was  believed  valid  for  two  reasons. 

First,  because  the  study  was  to  investigate  the  case  where 
many  PRNs  occupy  a  single  pixel,  the  averaging  techniques 
would  account  for  the  partial  portions  of  MTA  values  with¬ 
out  causing  scintillation.  Second,  for  the  fields  of  view 
under  investigation,  the  pixel  texture  areas  were  very 
long  narrow  figures  closely  approximating  rectangles  ex¬ 
cept  at  the  boundaries.  The  boundaries  were  edge  smoothed 
to  account  for  these  non-linearities. 


2.1.6  AVERAGE  TONAL  VALUE  ESTIMATION 

For  the  case  where  many  PRN  values  were  averaged  to  form 
a  pixel  tone,  an  algorithm  was  developed  to  provide  a  com¬ 
binatorial  process  that  could  be  implementable  in  suitable 
hardware  for  real-time  operations.  This  was  accomplished 
using  linear  estimation  theory.  Higher  order  estimators 
were  not  considered  as  their  hardware  counterpart  may  not 
be  suitable  for  real-time  operations. 

Let  Si j (n)  denote  a  "tonal  group"  or  matrix  of  n  x  n  MTA 
elements,  at  a  particular  location  (i,j)  in  the  PRN  data 
base.  Base  element  pij  is  defined  as  the  element  in  the 
upper  left  corner  of  Sj[ .  (n)  .  The  average  tone  of  Sij  (n) 
is  denoted  by  Aij(n).  3 

Since  direct  evaluation  of  Ai  j (n)  may  be  time  consuming 
for  large  values  of  n,  the  approach  here  is  to  seek 
approximations  to  Aij(n)  of  the  form 

Ai  j (n)  =  a  Pi  +  b 

where  a  and  b  are  the  correlation  coefficients  operating 
on  a  single  element  pi  of  the  Sij (n)  tonal  group.  The 
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element  pi  denotes  a  unique  element  of  the  tonal  group 
set  {pi-j/ .  These  linear  estimators  are  derived  in  the 
sense  or  classical  least  squares  estimators.  It  will  be 
demonstrated  that  these  estimators  do  exist  and  have 
several  remarkable  properties  within  the  set  of  all  maximal 
length  polynomials  of  degree  m. 

Since  S^.(n)  is  symmetric,  there  exist  only  K  =  2n  -  1 
unique  average  tonal  values  for  each  choice  of  n.  This 
reduced  set  of  ^ (n)  will  be  denoted  by 

Ai(n) 

The  classical  least  squares  estimator  will  be  derived  for 
n  =  2. 

Since  n  =  2, 


Ai(2) 


Pi  +  2Pi+l  +  Pi+2 

4 


The  difference  between  the  actual  average  tonal  value  and 
the  estimated  tonal  value  is  defined  as 


eik  -  V2>  “  <ak  Pi+k  +  bk> 

The  sum  of  the  squares  of  the  differences  is  expressed  as 

sk  - 1  c2ik  - 1  ‘V2>  -  ak  -  bk>2 


The  more  general  form  p^+>  is  included  here  since  the 
correct  element  of  the  code  set  {p^}  for  estimation  may 
not  be  the  base  element  of  S^-(n).  The  complexity  intro¬ 
duced  by  the  form  pi+.  now  leads  to  performing  k  regressions 
instead  of  just  one  regression.  Also  the  value  of  k  must 
be  determined  in  addition  to  the  constants  a^  and  b^. 


The  conditions  that  6^  be  a  minimum  are 

9  a. 

*  \  Ai  '  ak  2  Pi+k  -  "  bk  '  0 

and 

J  Ai  pi+k  "  ak  l  pi+k  _  bk  J  pi+k  *  0 
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This  system  of  equations  may  be  written  as 


?  Pi+k  nl 

~v 

r  ?  a.  i 

1 

i  i 

„  2 

m 

Pi+t  ?  Pi+k 

a. 

E  .  - 

Ll  l 

_  k  _ 

L  i  Pi+k  Ai. 

Let  C  be  the 

The  solution 

a,  =  (I 
K  i 

bk  =  (E 
K  i 

where  det  C  = 


coefficient  matrix  of 


set  is  obtained  using  Cramer's  Rule: 


Al  '  i  Pi+k  ‘  "  i  Pi«  ‘  Ai,/det  ° 

Pi+k  '  l  Pi+k  •  Ai  -  l  P\+k  *  £  Ai»/det  C 

(£  Pi+k*  -  n  E  p  i+k 


For  all  cyclic,  maximal  length  polynomials  of  degree  m, 
det  C  is  invariant  since  the  summations  include  all  members 
of  the  code  set  (pi). 

Since  the  estimator  for  A^d)  is 
A^l)  *  a  •  Pi  +  b 

where  a  =  1 
b  =  0 

the  correct  estimator  for  Adn),  n  >  1  will  be  found  from 
the  set  tafcr bjJ  for  that  value  of  k  for  which  a*  is  closest 
to  1  and  b£  is  closest  to  zero. 

An  examination  of  all  values  of  a^,  b^  for  all  eighth-degree 
polynomials  indicates  that  the  correct  value  of  k  is  that 
which  causes  Pi+fc  to  be  the  element  of  the  secondary  diag¬ 
onal  of  Sij(n).  This  is  consistent  with  intuitive  expec¬ 
tation  since  the  contribution  of  the  secondary  diagonal  is 
weighted  most  heavily. 
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There  are  eight  maximal  length  polynomials  of  degree  8. 

For  each  polynomial  a  search  for  a  best  estimator  yielded 
a  and  b  coefficients  that  were  identical  for  each  tonal 
order  n  of  S^j (n) .  This  means  that  the  estimation  pro¬ 
cedure  developed  for  this  simulation  is  independent  of  the 
PRN  code  or  textural  definition  of  the  data  base.  Other 
codes  providing  different  textural  detail  could  be  run 
using  the  same  simulation  software. 


2.2  DATA  I/O  AND  SIMULATION  CONTROL 

The  simulator  was  designed  to  accept  input  data  from  sev¬ 
eral  input  media:  cards,  disk  or  an  interactive  display 
terminal.  The  easiest  method  (and  most  popular)  was  found 
to  be  input  via  the  CRT  terminal  using  a  video  interactive 
menu  with  system  prompts.  The  program  operator  need  only 
respond  to  questions  displayed  on  the  CRT  terminal. 

The  data  input  routine  enables  the  operator  to  have  complete 
control  over  the  simulation.  The  data  were  divided  into 
four  main  areas:  data  base  parameters,  simulation  para¬ 
meters,  video  display  parameters,  and  aircraft  parameters. 

The  significant  parameters  which  are  under  user  control 
include: 

•  DATA  BASE 

Polynomial  Generator  -  PRN  size,  and  configuration 

Shading  -  Intensities  of  sky  and  ground 

Number  of  intensities  within  data  base 
Center  or  average  intensity  of  data  base 

Terrain  Description  -  Terrain  Size 

Terrain  Resolution 
Subterrain  Structure 


•  DISPLAY 

Size  of  sensor  display  screen 
Raster  display  resolution 

•  SIMULATOR 

Time  Scaling  -  Refresh  rate  of  screen 
Length  of  flight 

Scale  Factors  -  Integer  mathematical  routines 
Internal  Smoothing 
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Printout  Controls  -  Control  over  which  variables 

are  to  be  printed 

Rate  at  which  reports  are 
generated 


e  AIRCRAFT 

Position,  Velocity,  Attitude  (roll,  pitch,  yaw) 

These  variables  provide  the  user  with  a  method  for  describ¬ 
ing  the  total  working  environment  of  the  software  and  simu¬ 
lation  scenario. 

The  simulator  outputs  to  two  primary  units;  video  data  are 
put  onto  nine-track  magnetic  tape,  and  the  reports  are 
stored  on  an  online  disk.  The  tape  data  are  in  a  format 
compatible  with  the  video  line  assembler  of  the  ACTES  facil¬ 
ity.  Reports  are  kept  on  disk  primarily  to  increase  simu¬ 
lator  efficiency  and  to  avoid  tie-ups  on  the  line  printer. 


2.3  SPATIAL  TRANSFORMATIONS 

As  illustrated  in  Figure  4,  the  display  window  coordinate 
system  is  defined  by  the  right-hand  coordinate  system,  u,  v, 
and  w.  The  u  axis  is  perpendicular  to  and  intersects  the 
display  plane  at  its  center.  The  v  axis  is  in  the  direction 
of  scanning  a  display  line,  and  the  w  axis  is  in  the  dir¬ 
ection  of  successive  scan  lines  for  a  TV  raster  type  display 
The  origin  of  the  u,  v,  w  coordinate  system  represents  the 
location  of  the  aircraft  sensor  display  focus  or  view¬ 
point. 

The  fixed  world  coordinate  system  is  designated  by  the  right 
hand  coordinate  axes  x,y,z.  it  is  necessary  in  the  simu¬ 
lation  to  represent  a  point  (x,y,z)  in  the  fixed  world  sys¬ 
tem  from  the  viewpoint  of  the  aircraft?  that  is,  to  deter¬ 
mine  the  u,  v,  w  coordinates  of  the  point  (x,y,z).  To 
accomplish  this,  the  origin  of  the  x,y,z  system,  (Xq,y0,z0), 
is  first  translated  into  (x^,  y^,  z£,  which  is  the  dis¬ 
placement  of  (x0,y0,Zp)  in  respective  axis  directions  to 
the  location  of  the  aircraft  sensor  display  focus.  Then, 
for  any  point  (x,y,z)  a  point  (x^y'jZ')  is  found  relative 
to  (Xq,  y£,  Zq) • 
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Display  Window  and  Fixed  World  Coordinate  Systems 


Lastly,  to  account  for  rotation  of  the  aircraft  coordinate 
system  relative  to  the  fixed-world  system,  the  matrix 


x' 

y ' 
2* 


is  pre-multiplied  by  a  coordinate  transformation  matrix: 
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where  the  components  of  A  are  the  direction  cosine  relation¬ 
ships  with  respect  to  three  independent  angles:  heading, 
pitch,  or  flightpath  angle,  and  roll. 

The  point  u,v,w  is  the  three-dimensional  image  space  value 
which  is  to  be  clipped,  if  necessary,  and  mapped  into  the 
two-dimensional  sensor  image  plane. 

Point  ( Xq , y£ , Zq )  represents  the  location  of  the  viewpoint 
with  respect  to  the  origin  (Xg,y0, zQ)  bf  the  fixed-world 
system.  Aircraft  position  will  vary  from  TV  frame  to  TV 
frame;  thus,  for  each  frame,  a  new  position  (x£,y£,z")  is 
calculated  from  the  previous  position  (x£,y£,z£)  as  follows: 

x£  «  x£  +  vx  dt 

yo  *  yo  +  vy  dt 

Zo  “  zo  +  vz  dt 

where  vx,  vy,  vz  are  aircraft  velocity  vectors  in  the  x,y 
and  z  directions,  and  dt  is  the  time  interval  between  TV 
frame  calculations. 
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2.4  CLIPPING 


In  the  clipping  algorithm,  the  terrain  is  represented  by 
an  ordered  sequence  of  points;  consecutive  pairs  of  these 
points  comprise  terrain  edges.  Each  time  the  terrain  is 
clipped  against  a  particular  clipping  plane,  a  new  sequence 
of  points  results. 

The  mathematics  of  computing  a  point  of  intersection,  P:, 
between  a  clipping  plane  and  a  terrain  edge  (represented 
in  the  three-dimensional  coordinate  system  of  the  aircraft) 
is  as  follows: 

Pi  *  (ui'  V  wi> 

P^  *  (u. ,  v^,  w^)  =  first  endpoint  of  terrain  edge 
P2  *  ^u2'  v2 *  w2^  =  secon<*  endpoint  of  terrain  edge 

|PlPi|  =  distance  from  P^  to  P^ 

|P1P2|  *  distance  from  to  P2 

u±  *  ux  +  r(u2  -  ux) 

vi  ~  V1  +  r(v2  -  vx) 
w±  »  wx  +  r(w2  -  wx) 

lplpi|  I PiPi I 

where  0  <  r  <1,  and  r  =  *  p— 1  *  si  pi 

-  -  1*1*21  I *1*2 | 

P£  and  P2  are  illustrated  in  Figure  5a  and  defined  below: 

P1  *  <ui*  *  projection  of  Pj^  onto  u-v  plane 
P2  *  (u2*  v£)  *  projection  of  P2  onto  u-v  plane 

For  the  right  clipping  plane  illustrated  in  Figure  5b,  r  is 
computed  using  basic  geometric  relationships  to  be  the 
following: 


lvl  “  ui  tan&* 

r  “  lvi"ui  tanp|+|v£-u£  tane| 


idler e  5  *  1/2  horizontal  FOV 


-w 


Right 


Figure  5a.  Projection  of 

P1'P2  onto  u“v  plane 


Figure  5b.  View  onto  u-v 
Plane 


Similarly,  for  the  left  clipping  plane  r  can  be  determined 
from: 

|v^  +  u^  tan  & | 

r  *  lvi  +  ui  tan  P I  + lv2  +  u2  tan  P I 


For  top  and  bottom  clipping  planes,  respectively, 

|wj  +  uj  tan  a | 

r  =  |w|  +  u£  tan  of  +  w2  +  u2  tan  “I  *t0p  plaRe| 

|w^  -  u^  tan  a  | 

r  *  T*[~-~u;~tan  af  T  W2  -  ~u]  tan  d\  (bottoIn  plane) 

where 

a  *  1/2  vertical  FOV, 

Pi  *  (ui»wp  “  projection  of  P^  onto  u-w  plane, 
and  P£  ■  (u^wp  ■  projection  of  P2  onto  u-w  plane 

As  a  result  of  clipping,  the  vertices  which  define  the 
terrain  in  the  three-dimensional  viewpoint  coordinate  sys¬ 
tem  all  lie  within  the  field  of  view.  Thus,  they  can  be 
projected  onto  the  plane  of  the  screen  without  exceeding 
the  limits  of  the  screen  itself. 
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Scan  lines  are  numbered  from  -239  at  the  top  of  the  screen 
to  +240  at  the  bottom,  and  correspond  to  the  I  coordinate 
of  the  screen  points  (see  Figure  6) .  Horizontal  elements 
are  numbered  from  -320  at  the  left  of  the  screen  to  +320 
at  the  right  and  correspond  to  the  J  coordinate  of  the 
screen  points. 


J*  -320  J*  320 

1=  -240  I*  -240 


J=  -320  J*  320 

1=  240  I*  240 


Figure  6.  Display  Plane  Viewed  from  Focus  or 
Observer's  Position 


The  equations  which  project  a  point  (u'^v'^w')  onto  the 
screen  are  as  follows: 


A  W 

d  *  — 
u 


d  •  - 
u 


| horizontal  field  of  viewj 


where  d  ■  320  •  cot  y - j- 

■  distance  from  viewer  to  screen  in  screen  units 
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After  determining  the  screen  coordinates  of  terrain  verti¬ 
ces,  the  slope  and  J-intercept  of  each  screen  edge  can  be 
found.  This  information  is  needed  in  determining  where 
edges  are  crossed  by  scan  lines. 

^2  ”  ^1 
Slope  =  f — Z-T- 

I2  il 

Intercept  -  J2  ~  (sloPe)  '  X2 
where  and  (^f  J2)  are  endpoints  of  an  edge. 


2.5  PIXEL  TONE  CALCULATION 

The  primary  purpose  of  this  function  is  to  average  the 
appropriate  MTA-PRN  tones  which  map  into  a  single  picture 
element  or  pixel. 

Tonal  calculation  is  performed  when  part  (or  all)  of  a  TV 
line  falls  within  the  terrain  data  base;  it  is  performed 
on  a  line-by-line  basis.  The  variables  which  the  function 
requires  for  its  calculations  are  set  up.  These  include 
the  screen  and  ground  coordinates  of  the  first  pixel  to 
fall  within  the  terrain  boundaries,  the  location  and  size 
of  the  terrain  area  contained  in  the  pixel,  the  screen 
location  of  the  last  pixel  to  be  averaged,  and  edge  smooth¬ 
ing  constants. 

Once  these  variables  have  been  constructed,  an  algorithm 
averages  the  various  PRN  tones  which  fall  into  one  pixel. 
After  a  tone  for  an  individual  pixel  has  been  computed, 
it  is  scaled  and  stored  in  the  TV  data  line  in  its  appro¬ 
priate  location.  The  process  is  repeated  for  each  pixel 
on  the  TV  line. 

When  one  TV  line  has  been  computed,  edge  smoothing  of  the 
terrain  is  performed  and  the  final  tones  stored  for  print¬ 
ing  if  the  print  control  parameters  are  satisfied. 


2.5.1  PIXEL-TERRAIN  GEOMETRY 

Figure  7  shows  the  coordinate  systems  used,  the  fixed- 
world  coordinate  system  (x,y,z)  which  defines  the  terrain 
location,  and  the  viewpoint  coordinate  system  (u,v,w)  which 
defines  the  viewpoint  locations.  As  shown  in  Figure  7, 
the  (u,v,w)  system  is  a  translation  where  u  is  parallel 
to  x,  v  is  parallel  to  y,  and  w  is  parallel  to  z.  The  ro¬ 
tated  screen  coordinate  system  specified  as  (u',v',w') 
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ted 


Figure  7.  Simulator  Coordinate  System  Definitions 


where  the  u*  axis  i£  perpendicular  to,  and  passes  through, 
the  geometric  center  of  the  screen,  v'  is  parallel  to  the 
scan  lines,  w'  is  parallel  to  the  screen  vertical  axis. 


Trans formation  between  the  rotated  screen  coordinate  system 
and  the  unrotated  sensor  coordinate  systems  is  accomplished 
using  the  equation: 


V 

u 

V' 

= 

A 

V 

w’ 

m  - 

w 

where  the  components  of  [A]  are  the  directional  cosine 
relationships  with  respect  to  the  three  independent  angles 
of  heading,  pitch,  and  roll. 

Mapping  of  pixel  locations  to  terrain  coordinates  is  accom¬ 
plished  via  vectorially  projecting  image  plane  coordinates 
onto  the  flat  earth  using  the  equations 


xp  * 


•  Vup 


y  =  .  v  Ar 

P  P  P  P 

where  z  is  the  aircraft  altitude,  and  u  ,  v  and  w„  are 
P  P  P  P 

the  coordinates  of  a  corner  of  a  pth  pixel  in  the  (u,v,w) 

system. 


2.5.2  TONE  COMPUTATION 

Tonal  groups  have  been  developed  which  consist  of  N  x  N 
MTA  elements.  The  average  value  of  the  tonal  group  is 
estimated  using  a  linear  estimator,  which  reduced  compu¬ 
tation  time  considerably  at  a  minimum  loss  of  accuracy. 

To  accommodate  image  plane  rotation,  tonal  groups  will  be 
shifted  to  achieve  the  optimum  MTA  to  pixel  mapping. 

Figure  8  illustrates  the  two  basic  tonal  group  geometries. 
Figure  8a  shows  tonal  groups  which  are  stacked  vertically 
and  permitted  to  shift  horizontally  ("Free-Y"  tonal  groups) . 
Figure  8b  shows  groups  which  are  stacked  horizontally  and 
shifted  vertically  ("Free-X"  tonal  groups) .  This  descrip¬ 
tion  of  the  algorithm  will  examine  Free-Y  tonal  group 
development.  Free-X  groups  developed  utilize  the  same 
algorithm. 

The  size  and  number  of  tonal  groups  are  determined  using 
the  relationships: 
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Pixel  Projection 


Fixed  World  Coordinate 
System 


Figure  8a 


8a.  Tonal  groups  which  are  shifted  horizontally  (Free-Y) 
8b.  Tonal  groups  which  are  shifted  vertically  (Free-X) 


Figure  8.  Tonal  Group  Mapping 
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and 


AY, 


TG 


^TG 


AY 


PIX 


*PIX 


AX, 


TG 


^Stg 


AY 


PIX 


*PIX 


where 


AYtg;  AXtg  are  the  linear  distances  covered  by  the 

tonal  groups  which  will  map  into  the  pixel 
projection. 

AYpix  AXpix  are  the  linear  distances  covered  by  the 
'  pixel  projections. 


*tg 

*PIX 


is  the  area  of  the  tonal  groups . 
is  the  area  of  the  pixel. 


Pixel  projections  consist  of  M  tonal  groups  comprised  of 
N  x  N  MTA  elements;  therefore. 


aytg 

=  N  • 

LMTA 

axtg 

=  M  • 

N  •  LMTA 

^Stg 

=  M  • 

N2  •  LMTA' 

where  LMTA  is  the  length  of  one  side  of  an  MTA  element. 

From  these  relationships,  the  size  (N)  and  number  (M) 
tonal  groups  are  found  to  be 


N  -  Apix  /  (LMTA  •  AXpIX) 

M  -  ApIJ{  /  (LMTA  •  A Ypix  *  N) 
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Vertical  placement  of  tonal  groups  is  determined  as 
follows.  The  total  height  of  the  comhined  tonal  groups 
is  computed.  This  value  is  subtracted  from  AX.  This 
value  represents  the  difference  in  the  vertical  distance 
covered  by  the  pixel  projection  and  the  tonal  groups. 

This  difference  is  divided  evenly  between  the  top  and  bot¬ 
tom  of  the  pixel  projection,  yielding  the  vertical  position 
of  the  tonal  group  stack.  This  fixes  the  vertical  (X)  po¬ 
sition  of  every  tonal  group  as  each  must  lie  directly  on 
top  of  the  one  below  it. 

Each  individual  tonal  group  is  then  permitted  to  float  in 
the  Y  direction  (thus  the  name  "Free-Y*')  to  best  align  it 
within  the  pixel  projection.  The  Y  coordinate  is  computed 
by  substituting  the  X  coordinate  of  the  center  of  the  tonal 
group  into  the  equation  of  the  pixel  centerline. 


2.5.3  EDGE  SMOOTHING 

Edge  smoothing  is  called  on  a  line-by-line  basis  to  compute 
tones  of  pixels  which  contain  both  terrain  and  ground.  All 
computations  are  done  using  image  plane  coordinates  rather 
than  the  fixed-world  (terrain)  coordinates.  The  equations 
of  the  terrain  edges  on  the  screen  are  set  up  for  each 
frame  prior  to  any  tone  computations  and  stored  for  their 
use  here. 

The  tones  of  these  pixels  are  computed  using  the  tone  x 
area  relationship; 


Tone  *  ygnd  •  IGND  +  yter  •  TER 


idle  re 
Tone 

YGND 

yter 

IGND 

TER 


new,  smoothed  tone  for 
fraction  of  pixel  that 
fraction  of  pixel  that 
tone  of  the  ground 
tone  of  the  terrain  in 


a  pixel 

contains  ground 
contains  terrain 

that  pixel 


The  fractions  Y^p  and  YTgp  are  computed  using  geometric 
and  trigometric  relationships  of  a  line  (the  terrain  edge) 
cutting  a  square  (the  pixel) .  Floating  point  math  is  em¬ 
ployed  for  maximum  accuracy. 
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2.6  ADDITIONAL  CONSIDERATIONS 

Additional  functions  were  established  to  enhance  simulator 
operation.  They  are  the  use  of  integer  mathematics  for 
speed  enhancements  and  gray  scale  attenuation  to  provide 
control  of  the  texture  contrast  levels.  The  following 
paragraphs  describe  these  functions. 

An  analysis  of  the  simulator  timing  model  showed  that 
floating  point  instructions  in  various  locations  of  the 
software  were  increasing  the  run  time  of  the  simulator 
considerably.  Integer  mathematics  were  then  instituted 
as  follows. 

The  floating  point  numbers  were  pre-scaled  using  a  scale 
factor  (SFACT)  in  which: 

SFACT  (Sacle  Factor)  =  2RASSUB  =  SHIFT  (1,  RASSUB) 


where 


RASSUB  is  am  input  parameter,  and 

SHIFT  is  a  function  by  which  the  first  argument  is 
shifted  to  the  left  a  number  of  bits  equal  to 
the  second  argument. 

More  precisely,  A  =  SHIFT  (M,N)  -  M  *  2N  or 

SFACT  =  SHIFT  (1, RASSUB)  =  1  •  2RASSUB 

After  scaling,  mathematics  is  performed  in  an  integer  mode, 
then  numbers  are  restored  to  their  original  dimension. 
Restoring  is  accomplished  {in  the  general  case)  by: 

A  =  SHIFT  (M,-N)  =  M  •  2~N 

A  negative  value  for  the  second  argument  indicates 
a  shift  to  the  right. 

Gray  scale  reduction  is  an  enhancement  developed  to  pro¬ 
vide  user  control  to  attenuate  the  texture  contrast  viewed 
in  any  given  scenario.  This  control  is  also  needed  to 
permit  various  polynomial  degrees  to  be  used  without  major 
changes  to  the  software.  The  ACTES  display  is  capable  of 
256  gray  scales,  and  an  eighth-degree  polynomial  has  255 
individual  tones  so  that,  if  desired,  all  tones  could  be 
produced  on  the  display. 
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The  algorithm  associates  actual  tones  into  groups,  where 
the  number  of  groups  is  an  input  parameter.  The  tones  are 
linearly  scaled  into  these  groups  by  simple  division. 

Since  this  process  had  to  be  accomplished  for  every  pixel 
tone,  binary  division  (shifting)  was  implemented  on  the 
algorithm.  The  tones  were  shifted,  leaving  enough  high 
order  bits  to  associate  them  with  the  correct  tonal  group. 
This  allows  the  full  range  of  PRN  tones  to  be  attenuated 
about  a  central  value  in  increments  of  2,  4,  8,  16,  etc. 
Other  PRN  generations  having  more  or  fewer  tones  than  the 
display  can  be  handled  by  this  algorithm. 


2.7  VIDEO  LINE  ASSEMBLY 

Data  for  a  TV  frame  display  are  calculated  on  a  line-by¬ 
line  basis.  For  TV  lines  which  correspond  to  sky  or  back¬ 
ground  only,  two  buffers  containing  the  appropriate  tones 
are  initialized  at  the  start  of  the  simulation.  For  each 
TV  line  which  crosses  the  terrain  area,  pixel  intensity 
calculations  are  performed  and  overlay  the  640-element 
buffer  preset  to  the  ground  tone. 

Every  line  buffer  is  prefixed  by  a  frame  number  and  line 
number,  and  is  passed  to  a  software  video  interface  pro¬ 
vided  by  AFHRL  personnel.  The  interface  routine  converts 
data  words  received  from  the  Sigma  5  to  a  format  acceptable 
to  the  ACTES  facility  for  offline  video  display  generation 
and  buffer  storage. 
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3.0  RESULTS 


The  ability  to  produce  moving  pictures  with  realistic 
textural  cues  is  an  important  requirement  of  training 
programs  and  target  discrimination  applications.  The 
objective  of  this  study  was  to  investigate  a  texture- 
producing  technique  using  PRN  code  states  to  represent 
texture  tones  over  a  flat  area  of  terrain.  A  software 
model  was  constructed  to  simulate  hardware  algorithms  for 
computing  pixel  tones. 

The  software  simulation  has  been  demonstrated  to  be  an 
effective  tool  in  generating  a  variety  of  controllable 
textures  for  CIG,  which  are  spatially  correlated  to  insure 
realistic  displays.  Moreover,  the  operating  simulation 
provides  flexibility  and  interactive  control  to  a  user  to 
explore  and  evaluate  the  effects  of  modifying  the  texture 
data  base  definition  and  gray  scale  attenuation  of  a  dis¬ 
played  texture. 

Using  the  software  simulator,  scenarios  were  produced  to 
simulate  views  from  a  sensor  moving  over  surfaces  of  vary¬ 
ing  textures.  It  was  shown  that  as  the  sensor  moves  closer 
to  the  surface,  increased  texture  detail  and  contrast  are 
observed.  Dynamic  imagery  demonstrated  smooth  movement 
of  the  viewpoint  over  the  terrain. 

The  significance  of  the  simulation  is  that  a  technique  for 
producing  spatially  correlatable  textures  has  been  vali¬ 
dated.  Moreover,  the  algorithms  developed  can  be  implemen¬ 
ted  in  hardware  for  real-time  applications. 


3.1  EVALUATION  OF  PRN  CODE  TECHNIQUES 

The  principal  effort  of  the  study  was  to  develop  and  in¬ 
vestigate  the  combinatorial  algorithms  in  which  large  num¬ 
bers  of  MTA-PRN  tones  are  averaged  to  form  a  pixel  intensity. 
An  alternate  means  to  summing  and  averaging  each  individual 
element  was  analyzed  and  subsequently  implemented  in  soft¬ 
ware.  A  stepped  estimator  was  selected  for  investigation 
as  it  represented  a  simple  algorithm  with  minimum  compu¬ 
tation.  A  stepped  estimator  is  a  linear  estimator  in 
which  the  regression  coefficients  are  fixed  (i.e.,  A**1,B=0)  . 

The  basis  for  selecting  a  candidate  technique  was  that  it 
should  be  amenable  to  hardware  suitable  for  real-time 
operations.  Linear  and  stepped  estimators  are  of  this 
form.  The  linear  estimator  has  the  form 


37 


+  B 


T.  =  A  P 


i+k 


where  T . , . represents  all  the  combinations  within 

the  sequence 

p 

i+k  represents  a  PRN  tone  within  the  sequence 

where  k  is  the  displacement  from  the 
initial  starting  state 

A,B  are  regression  coefficients  of  the  best 

fit  estimator 

A  regression  analysis  was  used  over  the  complete  code  se¬ 
quence  to  determine  the  values  of  A  and  B  for  each  value 
of  k,  the  displacement  from  the  initial  starting  state. 

The  resulting  estimator  for  the  sets  of  values  (k,A,B) 
represents  the  best  linear  estimator.  It  enables  an  N  x  N 
texture  area  to  be  approximated  by  a  linear  combination  of 
a  single  PRN  value.  Statistics  for  the  estimator  were 
generated  to  compare  the  approximation  tone  to  the  actual 
tone  to  determine  the  best  values  of  A,B,  and  k.  Statis¬ 
tics  were  also  generated  for  the  stepped  estimator. 

Table  1  shows  the  results  of  the  linear  regression  analysis 
for  a  2  x  2  texture  area  for  one  of  the  eighth-degree 
primitive  polynomials.  The  k  value  is  given  in  column  DSP. 
The  values  in  column  LEVAR  are  a  measure  of  the  least  square 
error  between  approximating  tone  and  actual  tone  for  all 
2x2  sequences  using  the  kth  value.  A  and  B  are  the 
repression  coefficients  for  the  kth  value.  The  column  LEVAR 
shows  the  incoherent  properties  of  PRN  codes.  Beyond  the 
first  few  values  of  k,  the  error  reaches  a  maximum  and  re¬ 
mains  constant. 

Similar  regressions  were  run  for  other  N  x  N  texture  values 
to  develop  the  values  for  k,  A,  and  B.  It  was  shown  that 
the  off-diagonal  term,  k  =  N  -  1,  for  all  areas  of  the  form 
N  x  N  had  minimum  error.  These  values  were  programmed  in 
the  simulator. 

All  eighth-degree  maximal  length  polynomials  were  tested 
and  the  values  of  k.  A,  and  B  for  the  best  estimator  were 
identical.  Thus,  a  single  set  of  linear  coefficients  could 
be  applied  to  all  maximal-length  codes  of  a  given  degree. 
This  implies  that  the  hardware  algorithm  for  a  specific 
degree  PRN  polynomial  can  be  used  for  all  maximal  length 
polynomials  of  the  same  degree. 

The  linear  estimator  yields  a  hardware  implementation  of 
the  form  shown  in  Figure  9.  The  coefficient  KA  represents 
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TABLE  1.  SAMPLE  COEFFICIENT  SEARCH 


POLYNOMIAL  IS:  lOOlOlOtl 
SIZE  =  2X2 


DSP  LEVA* 

0  471982.1562 

1  1  3QQ.35-.  9,219 

2  471982,2500 

3  798392.4375 

1  877670.5625 

5  896325.4375 

'»  900399,5000 

I  901111.3750 

8  901107.9375 

9  900996. 7S00 

11  900950.7500 

12  900951.4375 

‘3 _ 9Q095Q. 750.0 

14  900951.2500 

15  900951.2500 

_ 16 _ SLQQ95-L.Q625 

1 7  900951.5000 

18  900951.3125 

20  900950.9375 

21  900951.1875 

32  900950 ,815.0. 

23  900951.2500 

24  900951.1875 

_ 25 _ 9M9iL.JLl5JI 

26  900951.6250 

27  900951.6250 

_29 _ 90  0_95 _L,  5  6.2  5. 

29  900951.2500 

JO  900951.3750 

31  900951.562 5 

32  900951.1250 

33  900951.8750 

34  _ 900951.3125 


A-COFFF 

0.557309 

0.747039 


0.557309 
0,272697 
JL.  1  303.39. 


0.059056 

0.023206 

0,004863 


0.005141 
0.010329 
0.011  BIX 


0.011811 

0.011811 

Q.Qliail. 


0,01181  1 
0,011811 
Q.  011  81 1 


0,011811 

0,011811 

(UOilHUL 


0,011811 
0.011811 
O.J 
0.011811 
0.01181 1 


0.01181 1 
0,011811 
SL.0.L1JLUL 
0.011811 
0,011811 
0.011811 


0.011811 
mum 


B-COFFF 

56.664421 

32.319051 


56.664421 
93.094780 
I  1 1.316628 


120.440872 
125.029686 
1 27.377487 


128,658081 
129,322098 
129.511 B1C 


129.511810 

129.511810 
_L2iL.BJJ.8Ul 


129.511810 
129.511810 
1 29.511  H1Q 


129.511810 
129.511810 
129.511 B10_ 


129.511810 

129.511810 
J2  9.5118 1  CL 

129.511810 
129,511810 

JL29^5J1JJXL 


129.511810 

129.511810 
-.129.  51181  _0_ 

129.511810 
129.511810 

_liSLa5.118l.flL 


129.511810 

129.511810 

.J2ft.MJ.ttUL 


attenuating  or  scaling  the  dynamic  range  of  the  PRN  values 
(i.e.,  255  shades  for  eighth-degree  polynomial),  and  Kq 
represents  a  D-C  offset  of  the  resulting  pixel  tone. 


3.2  PIXEL  TONE  ACCURACY 

An  analysis  was  developed  to  investigate  the  accuracy  of 
tonal  estimates  of  the  simulator.  Based  on  the  analysis, 
a  software  subsystem  was  implemented  which  produces  numeri¬ 
cal  results  of  the  estimations  used  in  the  simulator.  The 
software  package  computes  both  the  exact  tone  and  tone 
generated  by  the  simulator  algorithm  for  a  group  of  40  pix¬ 
els  located  in  the  center  of  a  TV  line.  These  two  tones 
are  compared  and  statistics  compiled  of  the  exact  versus 
simulated  tones  for  both  static  and  dynamic  scenarios. 

In  addition,  specific  analysis  software  has  been  implemen¬ 
ted  to  pinpoint  the  exact  degree  of  effects  of  algorithms. 


3.2.1  ANALYSIS  DESCRIPTION 

The  software  computes  the  actual  tones  for  a  group  of  40 
pixels  located  in  the  center  of  a  TV  line.  These  actual 
values  are  compared  to  those  tones  created  by  the  surface 
texturing  simulator.  Statistics  are  kept  on  the  difference 
between  the  actual  and  estimated  tones  of  both  static  and 
dynamic  scenarios  for  this  40-pixel  window. 

The  analysis  allowed  the  following  conditions  to  be  run: 

1.  Terrain  contrast  is  a  program  variable  to  show 
the  effects  of  higher  contrast  scenes. 

2.  The  analysis  generates  results  for  up  to  30 
consecutive  TV  lines,  each  40  pixels  wide,  to 
compare  line- to- line  variations  for  both  static 
and  dynamic  scenarios. 

3.  For  the  dynamic  scenarios,  the  program  collects 
the  results  of  consecutive  frames  for  the  same 
line  in  a  single  report. 

4.  All  of  the  8-bit  polynomials  used  in  the  simu¬ 
lation  could  be  run  to  develop  the  statistics 
for  the  different  PRN  generators. 
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Specific  analysis  are  also  performed  by  the  subsystem. 
These  analyses  include  pixel  projection  mapping  statis¬ 
tical  data  of  each  of  the  eight  PRN  sequences,  and  con¬ 
version  of  linear  estimator  coefficients  to  integer  values 


3.2.2  ANALYTICAL  RESULTS 

All  of  the  illustrations  and  reports  shown  has  been  genera¬ 
ted  using  the  same  parameters.  Aircraft  altitude  is  2400 
feet;  terrain  is  6  miles  down  range;  aircraft  velocity  is 
200  knots;  the  PRN  sequence  used  is  the  eighth-degree 
polynomial  x8  +  x5  +  +  x  +  1,  terrain  contrast  is  set 

to  32  shades,  and  the  TV  line  arbitrarially  selected  for 
reporting  statistics  is  line  364. 

The  analysis  included  investigations  of  the  following  top¬ 
ics: 

1.  Integer  Estimator  Coefficients 

2.  Static  Pixel  Statistics 

(a)  Theoretical  Pixel  Tone 

(b)  Static  Scene  Statistics,  Pixel  by  Pixel 

3.  Dynamic  Pixel  Statistics 

4.  Linear  Estimator  Local  Accuracies 


3. 2. 2.1  Integer  Estimator  Coefficients 

Within  the  surface  texturing  simulator,  the  coefficients 
of  the  linear  estimator  are  multiplied  by  2  7  or  128,  then 
converted  to  integer  values.  After  application  of  the 
linear  estimator  to  establish  a  tone  for  a  particular 
pixel,  the  tone  is  rescaled  by  multiplying  it  by  2~7. 

This  is  done  to  increase  the  speed  of  the  simulator. 

Comparison  tests  were  run,  correlating  the  estimated  tones 
of  tonal  groups  using  both  floating  point  coefficients 
and  scaled  integer  coefficients.  The  differences  between 
the  two  runs  were  negligible. 


3. 2. 2. 2  Pixel  Statistics 


One  of  the  capabilities  of  the  software  analysis  package 
is  to  compute  the  actual  tonal  value  of  a  given  set  of 
pixels  and  compare  this  tone  to  that  which  is  computed  by 
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the  surface  texturing  simulator.  The  actual  value  is 
computed  using  a  tone  x  area  algorithm  and  includes  all 
whole  »nd  any  fractional  parts  of  MTA  elements  which 
map  into  a  pixel.  This  value  is  compared  to  the  value 
which  is  computed  using  the  surface  texturing  simulator. 
Comparison  is  done  on  a  window  which  is  40  pixels  wide 
located  in  the  center  of  a  TV  line.  The  TV  line  is 
selectable  as  an  input  parameter.  The  program  will  com¬ 
pute  these  statistics  for  up  to  30  consecutive  TV  lines. 

Table  2  presents  the  results  of  this  analysis  phase. 

Columns  of  the  table  have  the  following  meanings: 

Pixel:  The  specific  pixel  of  the  40  the  row  of 

data  represents  (pixel  20  is  in  the  center 
of  the  screen) . 

256  Actual:  The  true,  unsealed  tonal  value  of  the  pixel 
element. 

32  Actual:  The  true  tonal  value  of  the  pixel  scaled  to 

32  shades. 

256  Est:  The  unsealed  estimated  tonal  value  of  the 

pixel  as  computed  by  the  sensor  simulator. 

32  Est:  The  estimated  tonal  value  scaled  to  32  shades. 

256  Del:  The  difference  between  the  actual  and  esti¬ 

mated  value  unsealed;  256  Actual  -  256  Est. 

32  Del:  The  difference  between  the  actual  and  esti¬ 

mated  value  scaled  to  32  shades;  32  Actual  - 
32  Est. 

256  Percent:  (256  Del/256)  *  100  (Percent  error  of  256 
shades) 

32  Percent:  Percent  error  in  32  shades  (32  Del/256)  *  100 


The  table  in  general  represents  a  slightly  worse  than 
average  case. 

The  columns  labeled  32  Del  and  32  Percent  are  the  most 
meaningful  as  they  represent  how  typical  scenarios  appear 
in  the  simulator.  The  least  accurate  pixel  on  the  page  is 
pixel  39  which  is  1.1%  or  three  shade?  off;  this  occurred 
in  less  than  4  percent  of  the  pixels  that  were  examined. 


TABLE  2.  SIMULATION  ESTIMATION  STATISTICS 


UNI 

NUM9ER 

144 

299 

12 

299 

92 

294  92  294 

92 

PIXEL 

ACTUAL 

ACTUAL 

1ST 

DEL  DEL  PERCENT 

PERCENT 

1 

149 

101 

114 

102 

12 

1  4.99 

0.99 

2 

192 

104 

141 

102 

11 

2  4.20 

0.79 

4 

190 

101 

144 

102 

9 

0  2.14 

0.00 

4 

149 

102 

192 

104 

-1  -1  -1.17 

•9.29 

9 

144 

101 

144 

102 

•5 

0  -1.99 

0.00 

9 

144 

101 

149 

102 

-4 

0  -1.94 

0.00 

7 

119 
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147 

109 

-4  -1  .9.92 

-0.29 

• 

144 
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101 

9 

2  9.92 

0.79 

• 
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122 

101 

7 

1  2.71 

0.99 

10 

119 

101 

122 

tot 

2 

0  0.74 

0.00 

11 

119 

102 

111 

101 

9 

1  1.17 

0.99 

12 

129 

101 

120 

100 

9 

1  9.92 

0.99 

11 

129 

100 

111 

94 

12 

1  4.49 

0.29 

14 

127 

100 

119 

94 

9 

1  9.12 

0.29 

19 

124 

100 

114 

94 

7 

1  2.79 

0.29 

14 

121 

100 

119 

99 

9 

1  2.94 

0.29 

17 

111 

99 

119 

99 

-9 

0  -1.95 

0.00 

If 

114 

99 

129 

100 

1  -9.52 

-0.24 

It 

117 

99 

112 

101 

*14  * 

2  -5.99 

-0.74 

20 

109 

99 

127 

100 

m  1  (| 

2  -7.09 

-0.74 

21 

107 

94 

129 

100 

*  f  4  * 

2  -9.29 

-0.74 

22 

107 

99 

111 

99 

-4 

0  -1.54 

0.00 

21 

111 

99 

119 

99 

•2 

0  -0.74 

0.00 

24 

109 

99 

112 

94 

-4 

1  -1.54 

-0.99 

29 

44 

97 

107 

99 

-4 

1  -2.92 

•0.99 

20 

92 

99 

111 

99 

**  1 9 

2  -7.42 

-0.74 

27 

100 

97 

109 

99 

-9 

1  -4.11 

-0.99 

21 

107 

99 

119 

to 

-11 

1  -4.20 

-0.99 

29 

111 

91 

129 

100 

*  1 9 

2  -4.99 

-0.79 

10 

104 

99 

120 

100 

•  1 2 

2  -4.40 

-0.74 

•  1 

114 

99 

119 

99 

i 

0  0.90 

0.00 

12 

117 

99 

119 

99 

2 

0  0.79 

0.00 

11 

120 

100 

129 

100 

■2 

0  -1.17 

0.00 

14 

129 

101 

124 

100 

9 

1  1.99 

0.99 

19 

119 

102 

120 

100 

19 

2  7.42 

0.79 

19 

140 

101 

124 

100 

19 

2  9.25 

0.79 

17 

111 

101 

122 

101 

2 

0  1.17 

0.00 

•9 

197 

102 

129 

101 

9 

1  9.19 

0.99 

19 

141 

102 

119 

99 

29 

4  4.77 

1.17 

40 

141 

102 

121 

100 

20 

2  7.91 

0.79 

3. 2. 2. 3  Dynamic  Pixel  Statistics 


Table  3  illustrates  a  dynamic  scenario  which  compares  actual 
to  estimated  tones  of  the  40-pixel  window.  The  numbers  are 
hexadecimal  representations  of  the  tonal  values.  Each  row 
pair  represents  one  frame  snapshot  of  the  40-pixel  window 
for  the  true  pixel  tones  (top  row)  and  the  estimated  tones 
(bottom  row) .  The  next  row  pair  shows  the  tonal  value  for 
the  next  successive  frame.  The  simulation  information  is 
the  same  for  this  set  of  numbers  as  those  for  Table  2.  The 
numbers  show  that  the  data  base  holds  consistent  over  a 
dynamic  sequence. 

The  dynamic  scenario  shows  the  same  statistics  as  the  static 
pixel  statistics.  There  is  a  maximum  of  ±2  shade  variations 
between  actual  or  exact  tone  and  simulated  tone. 


3.3  SIMULATION  RESULTS 

After  the  validity  of  program  operations  were  verified,  al¬ 
gorithms  were  optimized  for  improved  execution  time.  A 
timing  model  was  developed  for  the  simulator.  It  revealed 
that  subroutine  linkage  was  an  impedance  which  could  easily 
be  altered.  More  importantly,  it  showed  that  floating 
point  mathematics  instructions  in  several  of  the  averaging 
algorithms  were  consuming  the  most  time;  thus,  floating 
point  mathematics  functions  were  replaced  by  integer  mathe¬ 
matics  functions. 

These  two  enhancements  combined  to  give  a  savings  of  more 
than  50  percent  so  that  typical  run  times  were  in  the  order 
of  30  to  40  seconds  for  producing  a  TV  frame. 

A  six-second  scenario  was  produced  which  showed  the  terrain 
moving  smoothly  and  without  scintillation.  Of  note  was 
the  fact  that  the  simulation  produced  a  different  tone  per 
pixel  on  each  raster  line;  this  pointed  out  the  need  for  an 
extremely  high  quality  monitor  that  was  properly  calibrated 
as  to  TV  raster  and  interleaving  effects. 

For  speed  in  execution,  table  lookup  techniques  proved  to 
be  the  best  method  for  defining  the  basic  tones  of  the  PRN 
structure  and  starting  point  of  the  terrain  structure. 

This  technique  also  allowed  the  study  of  subterrain  struc¬ 
tures. 
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TABLE  3.  DYNAMIC  ACTUAL/ESTIMATED  COMPARISON 
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Texture  data  base  sizes  were  selected  such  that  the  terrain 
was  typically  limited  to  1.5  nautical  miles  in  length 
(dependent  on  MTA  sizes).  There  was  no  limitation  on 
terrain  width.  The  program  could  be  altered  to  expand  to 
larger  terrain  lengths  up  to  the  limits  of  available  user 
memory  on  the  operating  computer.  Typical  terrains  in¬ 
vestigated  were  between  600  feet  to  2  nautical  miles  wide 
and  1.5  miles  in  length. 


3. 3.1.1  Texture  Types 

Figure  10  is  a  set  of  still  photographs  that  show  different 
texture  patterns  using  the  same  PRN  polynomial.  The  photos 
were  taken  with  extremely  high  contrast  on  the  display  so 
that  they  could  be  reproduced;  the  actual  visual  display 
can  be  set  with  less  contrast,  which  is  more  pleasing  to 
the  viewer.  Some  effect  of  contrast  of  the  texture  will  be 
described  later  in  this  section. 

Table  4  defines  the  scenario-simulation  parameters  for 
Figure  10.  This  figure  demonstrates  the  versatility  of 
PRN  codes  as  a  source  of  texturing  data.  All  of  the  pic¬ 
tures  shown  in  the  figures  were  created  by  the  same  eighth- 
degree  polynomial  generator.  The  variation  that  is  shown 
depicts  various  forms  of  data  base  structure  and  subterrain 
structure.  Figure  10  (b-c)  show  mosaic-type  structures 
while  lOd  depicts  a  more  continuously  modulating  structure. 
It  is  believed  that  greater  variations  in  texture  type  can 
be  obtained  by  usage  of  different  and  larger  polynomial 
generators,  different  subterrain  structures  and  various 
shading  schemes.  A  detailed  investigation  of  the  various 
mosaic  and  continuous  terrain  types  was  beyond  the  scope 
of  this  study  but  is  recommended  for  future  study. 

Figure  10a  shows  the  basic  texture  created  by  an  8-bit  PRN 
polynomial  generator  with  the  stepped  algorithm  subterrain 
structure.  Here  the  subterrains  are  clearly  visible.  Re¬ 
view  of  the  photograph  will  show  the  basic  pseudo-random 
pattern.  Several  bands  are  wide,  while  others  are  very 
narrow.  Some  bands  change  gradually  from  dark  to  light 
(and  vice  versa)  while  others  change  abruptly.  In  this 
particular  photograph,  the  8 -bit  generator  does  not  repeat 
itself  in  the  horizontal  direction  but  does  repeat  itself 
several  times  in  the  vertical  direction. 

Review  of  the  photographs  on  the  actual  ACTES  display  shows 
the  repeating  pattern  in  the  terrain  background  but  with 
gradual  loss  of  foreground  detail.  The  photographs  in 
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Figure  10  (continued) 
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Figure  10  also  show  the  results  of  the  combinatorial 
algorithms  and  data  base  quantization.  Figure  10b  shows 
the  effects  of  data  base  quantization.  At  the  bottom  of 
the  terrain  one  MTA-PRN  tone  covers  more  than  one  pixel. 

At  the  top  of  the  terrain,  many  MTA-PRN  tones  are  averaged 
to  form  the  pixel  tone.  Depending  on  the  mission  flyby 
geometry,  the  minimum  texture  area  (MTA)  should  be  selected 
such  that  the  simulation  is  always  averaging  many  PRN  tones 
per  pixel  to  eliminate  a  stylized  display,  thereby  pro¬ 
viding  more  display  realism.  Intra-line  pixel  averaging 
could  also  be  implemented  to  smooth  this  effect  if  it 
proves  to  be  a  desirable  effect  for  terrain  cueing. 


3. 3.1.2  Texture  Contrast 

Figure  11  (a-d)  shows  a  terrain  at  close  range  (300  ft., 

2  nm)  with  different  gray  scale  levels  or  attenuation  of  the 
tonal  dynamic  range  within  the  terrain  structure.  Table  4 
describes  the  simulation  parameters.  The  original  255 
shades  generated  by  the  PRN  generator  are  scaled  down  to 
32,  16,  8  and  4  shades,  respectively,  for  the  four  photo¬ 
graphs.  It  is  easily  seen  how  the  texture  becomes  smoother 
or  undulating  as  the  number  of  shades  decreases.  Increas¬ 
ing  the  number  of  shades  results  in  a  more  detailed  terrain 
image  with  sharp  contrast  levels.  This  permits  greater 
flexibility  in  terrain  modeling  and  possibly  in  perception 
and  distance  cueing.  The  effect  is  that  as  the  number  of 
shades  increases,  the  terrain  appears  to  get  closer  to  the 
viewer  as  it  is  more  detailed,  regardless  of  the  fact  that 
the  images  do  not  increase  in  size. 


3. 3.1. 3  Texture  Detail 

One  of  the  most  important  concepts  shown  by  this  study  is 
the  effect  of  the  PRN  averaging  techniques  on  perception 
and  texture  cueing.  The  viewing  of  a  recognizable  "texture" 
pattern  which  increases  in  detail  as  the  range  from  viewer 
to  terrain  decreases  is  a  fundamental  to  flight  training 
and  target  discrimination  applications. 

Figure  12  (a  through  d)  shows  the  effect  on  textural  per¬ 
ception  as  the  terrain  approaches  the  aircraft  sensor. 

Figure  12a  shows  the  view  of  the  terrain  from  an  altitude 
of  1000  feet  and  range  of  5  nautical  miles.  In  Figure  12b, 
the  terrain  is  3  nautical  miles  in  front  of  the  aircraft  at 
the  same  altitude.  Note  that  in  the  texture  area  of 
Figure  12b,  the  texture  pattern  can  now  be  observed.  In 
Figure  12c,  the  aircraft  is  at  an  altitude  of  300  feet  and 
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32  Shades 


16  Shades 


Figure  11 


Altitude  =  1000  feet 
Figure  12 
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12c-  Range  =  2  miles 

Altitude  =  300  feet 


Figure  12  (continued) 
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a  range  of  2  nautical  miles  from  the  terrain.  From  this 
position  a  large  increase  in  texture  pattern  detail  can  be 
seen.  Table  5  summarizes  combinatorial  aspects  of  the 
simulation,  which  is  required  in  each  case  to  develop  the 
video  representation. 

The  relationship  of  texture  detail  to  distance  is  also 
illustrated  in  Figure  12d.  In  this  long  terrain,  the 
texture  pattern  is  visible  close  to  the  viewer  but  is 
lost  in  distance  with  only  the  general  texture  pattern 
apparent  to  the  viewer. 


TABLE  5.  AVERAGING  REQUIREMENTS  OF  PRN  ALGORITHM 


Figure 

Altitude/Range 

HTA's  Averaged 
(6*  Sq.  MTA's) 

■KTHH 

1000  ft  5  nn 

17 

86  X  3  -  258 

Hg  1 

1000  ft  3  nn 

184 

10 

31  x  2  -  62 

BB 

300  ft  2  nn 

Top  (Background) 

275 

7 

46  x  1  -  46 

Bottom  (Foreground) 

16.9 

1.7 

3x1-3 

3.3.2  FIELD  OF  VIEW  CONSIDERATIONS 

Figure  13  (a-c)  shows  the  effect  of  a  change  in  the  field 
of  view.  The  terrain  data  base  in  all  cases  is  exactly 
the  same.  It  is  positioned  the  same  distance  from  the  view 
point  in  all  instances.  Field  of  view  is  the  only  para¬ 
meter  which  is  varied.  Figure  13  (a  through  c)  represents 
fields  of  view  of  45  degrees,  20  degrees  and  10  degrees, 
respectively. 

Figure  13d  is  a  superposition  of  a  simulated  runway  (1.5 
nautical  miles  x  3000  feet)  over  a  textured  terrain.  The 
photograph  is  a  simple  illustration  of  a  textured  scene. 

The  runway  is  .25  nautical  mile  to  the  right  of  the  view¬ 
point. 

The  larger  fields  of  view  (greater  than  20*)  were  run  to 
demonstrate  the  applicability  of  the  technique  to  provide 
texturing  to  visual  CIG  displays.  Because  the  simulation 
algorithms  provided  the  same  types  of  presentations  as 
for  narrow  fields  of  view,  these  techniques  can  probably 
be  used  for  visual  presentations  in  wide  FOV  simulation 
applications. 
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13d.  Runway  in  a  Field 


Figure  13  (continued) 
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3.3.3  THREE-DIMENSIONAL  ORIENTATION 


The  simulation  was  enhanced  to  provide  six  degrees  of  free¬ 
dom  of  the  viewpoint  with  respect  to  the  terrain.  The 
terrain  edges  in  the  object  space  rotated  according  to  the 
angular  attitude  of  the  field  of  view.  The  texturing 
patterns  were  also  rotated  or  skewed  in  a  similar  manner. 


3. 4  ACTES  FACILITY  OPERATIONS 

The  simulation  was  delivered,  implemented  and  run  on  the 
Simulation  and  Training  Advanced  Research  System  (STARS) 
ACTES  Xerox  Sigma  5  computer  system  located  at  Wright- 
Patterson  AFB.  The  Sigma  5  computer  ran  the  actual  simu¬ 
lation  and  stored  the  video  data  on  magnetic  tape.  The 
offline  data  storage  retrieval  and  display  capabilities  of 
the  ACTES  facility  were  used  to  generate  and  display  the 
flyby  scenarios. 

The  generation  and  storage  of  video  data  from  the  Sigma  5 
onto  magnetic  tape  have  limitations.  Due  to  the  format 
required  for  the  ACTES  offline  capabilities,  approximately 
30  TV  frames  or  one  second  of  display  would  fit  on  a  tape. 
The  combined  lengthiness  of  the  computer  operations  (tape 
mounting,  rewind,  etc.)  impeded  the  generation  of  long 
simulations.  It  is  recommended  that  the  direct  Sigma  5  - 
ACTES  computer  interface  be  investigated  to  speed  up  the 
mechanical  process  of  generating  scenarios. 

Programs  were  initially  delivered  on  cards  and  stored  on 
disk.  Online  terminals  are  provided  for  source  file  up¬ 
dating,  program  compilations  and  simulation  data  input/ 
output  and  control  in  an  interactive  mode.  The  high  speed 
terminals  enabled  easy  access  to  source  files  and  made  the 
interactive  simulation  possible  and  simple  to  operate. 
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4.0  CONCLUSIONS  AND  RECOMMENDATIONS 
4.1  CONCLUSIONS 

This  study  showed  the  versatility  and  capabilities  of 
using  prn  coding  techniques  for  generating  texture  patterns 
in  CIG-type  simulations.  Variation  in  texture  detail 
was  demonstrated  and  showed  a  smooth  transition  from  fore¬ 
ground  detail  to  background  general  patterns.  A  number  of 
terrain  types,  including  mosaic- type  areas  and  continuous 
undulating  areas,  were  developed  showing  the  texturing 
capability  of  the  techniques.  Scenes  for  different  fields 
of  view  were  analyzed  and  indicated  that  the  technique  can 
be  extended  to  include  visual-type  displays  as  well  as  sen¬ 
sor-type  displays. 

The  techniques  implemented  in  this  non- real- time  software 
simulation  are  adaptable  to  real-time  hardware. 

During  playback  on  the  TV  monitor,  areas  of  scintillation- 
type  variations  were  noticed  in  both  the  terrain  texturing 
and  the  monotone  gray  areas.  Analytical  results  of  a  line 
of  TV  pixels  show  there  is  a  maximum  variation  of  i2  gray 
shades  within  the  textured  terrain  between  exact  and  simu¬ 
lated  tones.  It  is  believed  that  this  variation  would  be 
difficult  to  detect  within  the  texture.  Because  the  mono¬ 
tone  gray  area  had  similar  variations,  no  conclusions 
could  be  made  regarding  the  source  of  the  variations. 

Due  to  the  computation  time  required  to  create  a  scenario, 
long  scenarios  were  not  created  in  the  Sigma  5.  This  would 
have  assisted  in  the  determination  of  the  source  of  these 
variations.  It  is  recommended  that  a  partial  or  total 
hardware  implementation  of  the  tonal  calculation  be  made, 
thus  reducing  scene  computation  time  and  enabling  closer 
analysis  of  the  video  display  effects. 


4.2  SYSTEM  ENHANCEMENTS  AMD  RECOMMENDATIONS 

One  recommended  change  is  a  hardware  adaptation  to  the 
ACTES  system  to  perform  the  tonal  calculations  in  real  time. 
The  simulation  software  can  be  made  to  function  within  the 
ACTES. 

An  alternate  recommendation  is  a  partial  or  total  hardware 
implementation  in  a  high-speed  data  input /output  mode  of 
the  tonal  generation  presently  being  performed  in  software. 
This  implementation  will  increase  the  accuracy  and  speed  of 
the  tonal  calculations,  thereby  providing  the  capability 
to  produce  long  scenarios.  Long  scenarios  will  enable  a 
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visual  analysis  of  the  variations  presently  being  viewed 
on  the  TV  monitor.  This  hardware  could  become  a  model 
for  real-time  hardware  specifications. 

A  recommended  enhancement  to  the  system  is  the  ability  to 
provide  three-dimensional  texturing  in  the  object  space. 
Completion  of  this  enhancement  would  then  allow  the  use  of 
the  cyclic  coding  techniques  to  be  considered  for  inclusion 
in  existing  non-real-time  and  future  real-time  simulation 
systems. 

Further  investigation  of  cyclic  code  generators  should 
attempt  to  correlate  texture  patterns  to  realistic  terrain 
environments,  such  as  forests,  fields,  grass,  etc.  This 
effort  would  include  the  use  of  larger  code  generators, 
different  and  potentially  non-maximal-length  polynomials, 
subterrain  structures,  and  superimposed  multiple  code 
sequences. 
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